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This review concentrates on progress in lattice QCD during the last two years and, particularly, its im- 
pact on phenomenology. The two main technical developments have been successful implementations 
of lattice actions with exact chiral symmetry, and results from simulations with two light dynami- 
cal flavours which provide quantitative estimates of quenching effects for some quantities. Results 
are presented for the hadron spectrum, quark masses, heavy-quark decays and structure functions. 
Theoretical progress is encouraging renewed attempts to compute non-leptonic kaon decays. Although 
computing power continues to be a limitation, projects are underway to build multi-teraflops machines 
over the next three years, which will be around ten times more cost-effective than those of today. 



1 Introduction 

1.1 Overview 

Progress in lattice QCD tends to be in- 
cremental. Through improved formulations 
and increased computer power, we have been 
steadily gaining control over all the sys- 
tematic approximations inherent in numer- 
ical simulations. Rarely, progress is rev- 
olutionary. Happily, we are in the midst 
of such a major leap forward, through re- 
cent demonstrations that lattice formula- 
tions, which preserve exact chiral symmetry, 
work in practice. Combined with increasingly 
cost-effective computing technology, a period 
of accelerated progress, impacting directly on 
phcnomcnologically important QCD calcula- 
tions, can be foreseen over the next five years. 

In this review, I describe recent progress 
demonstrating that chirally symmetric for- 
mulations are feasible and results for phe- 
nomcnologically relevant quantities. I omit 
results for QCD thermodynamics and, with 
one exception, for non-QCD theories. Al- 
though there has been a lot of work on 
the challenging problem of the confinement 
mechanism, the understanding achieved so 
far is partial and the picture is still too con- 
fusing to do it justice in a review such as this. 
For results in these areas, other exploratory 
phcnomcnological applications and most of 



the technical details, I refer you to the lat- 
est in the annual series of lattice conference 
proceedings!!! 

1.2 Objectives of Lattice Field Theory 

The primary objective of lattice field theory 
is to determine the parameters of the Stan- 
dard Model and, thereby, to seek signals of 
new physics. Due to confinement, the quark 
sector is not directly accessible by experi- 
ment and numerical simulations of QCD are 
needed to provide the missing link. In prin- 
ciple, lattice QCD offers model-independent 
computations of hadronic masses and matrix 
elements. Ultimately, it should test QCD as 
the theory of strong interactions and provide 
an understanding of confinement. The name 
of the game is the control of systematics, par- 
ticularly to quantify dynamical quark effects 
and reliably simulate at, or extrapolate to the 
physical values of the light and heavy quark 
masses. 

The second objective is to determine the 
phase structure of hadronic matter. Both 
the location and the order of the line of 
phase transitions, separating the confined, 
hadronic phase from the deconfined, quark- 
gluon plasma phase, in the temperature (T), 
chemical potential (/z) plane, are sensitive 
to the flavour content. Since T c is close to 
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the strange quark mass, it is particularly im- 
portant to simulate the strange quark accu- 
rately, and this awaits simulations with more 
realistic dynamical light flavours. Also, the 
fi =/= plane is not accessible to simula- 
tions with current algorithms, because the 
action is complex and Monte Carlo impor- 
tance sampling fails. Significant progress 
has been made in QCD thermodynamics at 
/i = 0, and new approaches using anisotropic 
lattices should yield spectroscopic results to 
help better understand leptonic decays of 
vector mesons, strangeness production and 
J/tp suppression. 

Finally, looking beyond QCD, we aim 
to develop simulations into a general pur- 
pose non-perturbative tool. Recent progress 
in formulating lattice chiral symmetry has 
reawakened hopes of being able to simulate 
chiral and SUSY theories. 

2 Theoretical Progress 

2.1 Lattice Chiral Symmetry 

Major progress has been achieved over the 
past eight years following the rediscovery of 
the Ginsparg- Wilson (GW) relatione!. This 
states that if the lattice Dirac operator, D, is 
chosen to satisfy 



75D + Dj 5 = aDj 5 D, 



(1) 



where a is the lattice spacing, then the theory 
possesses an exact chiral symmetry. Such a 
formulation has the great virtue that it sep- 
arates the chiral and continuum limits. It 
also forbids 0(a) terms, so that the result- 
ing fermion actions are improved. The GW 
relation languished for many years, because 
there was no practical implementation. The 
breakthrough was the discovery of three con- 
structions: the overlapa, domain walla and 
perfect actional. As a consequence, lattice 
simulations with exact chiral symmetry are 
now a reality (see Kikukawa's talk! and other 
reviews!) . 



For vector theories like QCD, this means 
that we can maintain a global chiral symme- 
try at non-zero lattice spacing, so that simu- 
lations should be able to approach the phys- 
ical u and d quark masses in a controlled 
fashion. More importantly, the mixing of 
operators of different chiralities, which has 
plagued kaon mixing and decay calculations, 
is avoided. Formally, abelian and non-abelian 
chiral gauge theories have been constructed 
on the lattices, so the Standard Model can 
now be defined non-perturbatively. Some- 
what more speculatively, for SUSY theories 
without scalar fields, lattice chiral symmetry 
forbids relevant SUSY-violating terms in the 
action, and so offers the prospect of lattice 
simulations without fine tuninga. 

2.2 Overlap Quarks 

In the overlap formulation, the lattice Dirac 
operator has the form 



-Dqv(m) = 



1 + fJ, 



7 5 sgn(# w ), 
(2) 

where fi is the bare quark mass, and -f/w is 
the (hermitian) Wilson-Dirac matrix with a 
large negative mass. D ov (0) obeys the GW 
relation (Eq. (Q)). Numerically, the chal- 
lenge is to compute accurately the sign of the 
large sparse matrix. Hyj, eg using a rational 
approximation!!!]: 

N 



sgn(iJ w ) « ■ 



fur- (3) 

This approximation breaks down for small 
eigenvalues, so it is necessary to project out 
the lowest eigenvectors and treat their signs 
exactly. 

The resulting implementation has been 
tested for quenched QCDEl. The results in 
Figure [l] show that, as expected, there is no 
additive quark mass renormalisation, ie 

a = Z m V [1 + 0(a 2 )} . (4) 

The quenched light hadron spectrum has also 
been computed!!!!, albeit in a small volume, 
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Figure 1. The renormalised quark mass, obtained from the axial Ward identity, versus bare quark mpss 
(fj, = rriQa) in quenched QCD, computed using the overlap formulation at three different lattice spacingsEJ. 



and evidence is found that the mass ratios 
roughly scale at fixed quark mass, showing 
that it is now possible to do phenomenologi- 
cally interesting calculations in these new chi- 
rally symmetric formulations. 

Unfortunately, even for quenched QCD, 
the computational cost of simulating overlap 
quarks is comparable to that of dynamical 
simulations^, due to the extra work involved 
in computing the sign function (Eq. ([|)). Lit- 
tle is known yet about the cost of simulating 
dynamical quarks this way. 

The approximation of the sign function 
in Eq. ([|) may be replaced by local intera- 
tions amongst a set of 2N auxiliary fields: 

- N 1 

i) sgn(ffw) "0 ~ "0 ~ 

s=l c s + 

X! i(^ s Xs + Xsi>s) + Xs(c 2 s H w )Xs 

s 

+b s (xs<f>s + (f>sXs) ~ 4> s H W (l>s] ■ (5) 

This shows that the overlap formulation may 
be thought of as five-dimensional . in which 
the fifth dimension is like flavourEj. 

2.3 Domain Wall Quarks 

Here the fermions live in five dimensions and 
are coupled to a mass defect located on a 
four-dimensional hyperplane. On a finite lat- 
tice, with L s sites in the fifth dimension, the 
Dirac operator may be written 



-Ddw(m) = 




x 75 tanh ( -ylog Tj . (6) 

Here T is the transfer matrix in the fifth di- 
mension and (i is again the bare mass. The 
strong similarity with the overlap, Eq. (||), 
is evident. In fact, the two formulations are 
identical in the limits of L s — > oo and zero 
lattice spacing in the fifth direction. 

For finite L s , the chiral modes of opposite 
chirality are trapped on the four-dimensional 
domain walls at each end of the fifth dimen- 
sion. Chiral symmetry is broken, but the 
breaking is exponentially suppressed by the 
size of the fifth dimension. Sever al gro ups 
have tested this in quenched QCDtiH It 
is found that the pion mass does not always 
vanish with quark mass, in the limit L s — > oo, 
due to near unit eigenvalues of T, which al- 
low unsuppressed interactions between the 
LH and RH fermions. This is a strong- 
coupling effect, which goes away for weak 
enough coupling, or usine; a renormalisation- 
group improved actionEj. The problem may 
be controlled numerically by projecting out 
the low eigenvectors and taking their contri- 
bution with infinite L s t3. 

Thus, there is now a good understanding 
of how to achieve a close approximation to 
exact chiral symmetry in lattice QCD simu- 
lations. Although numerically relatively ex- 
pensive, it is early days and more efficient 
algorithms may yet be found. Already, the 
extra degree of control given by exact sym- 
metry probably outweighs the cost for matrix 
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Figure 2. String h*«aking in the four-dimensional 
SU(2) Higgs modeO. The results for the potential, 
V(r), versus separation, r, at two lattice spacings are 
in excellent agreement. 



element calculations, and we can anticipate 
rapid progress over the next few years. 

3 String Breaking 

For the remainder of this review, I will de- 
scribe results obtained using traditional for- 
mulations of lattice QCD, focusing on the ef- 
fects of dynamical quarks in order to try to 
quantify quenching errors for as many quan- 
tities as possible. 

Perhaps, the first effect of dynamical 
quarks we might hope to see is string break- 
ing. This flattening of the potential between 
two static charges, as they are separated, has 
been observed as level crossing in the confine- 
ment phase of the SU(2) Higgs modelO. At 
a particular separation, the string state be- 
comes degenerate with the state of two static- 
light "mesons" , as shown in Figure || 

Demonstrating string breaking in QCD 
at zero temperature is proving to be much 
more challenging than expected. As yet there 
is no completely convincing signal. This is 
because there is poor overlap between the 
string states used to compute the static quark 
potential and the broken-string state, com- 
prising two static-light mesons. Including the 
latter is computationally very costly, because 



it requires quark propagators at all sites. 

However, the mixing matrix element has 
been computed for two dynamical flavours 
and found to be non-zerdl3. Using only string 
states, there are hints of a flattening poten- 
tial, just as the signal becomes swamped by 
noiseEa, but a recent high-statistics calcula- 
tion provides pretty conclusive evidence that 
such attempts are doomedBa and the two- 
meson state must hii included (as was done 
for the Higgs modeltZl) . The situation is frus- 
trating, but it is only a matter of time before 
sufficient computing resources are brought to 
bear. 



4 Hadron Spectrum 

The important result that the quenched light 
hadron spectrum disagrees with experiment 
was finally established by the CP-PACS 
Collaboration^] in 1998 and announced at 
ICHEP98. This had proved difficult, requir- 
ing high statistics, because the deviation is 
less than 10%. This small deviation is good 
news for phenomcnological applications of 
lattice QCD, which still rely heavily on the 
quenched approximation. The main symp- 
tom of quenching is that it is not possible 
consistently to define the strange quark mass 
- the two spectra obtained from using the K 
and the <fi meson to determine the strange 
quark mass disagree. 

Since 1998, the focus has been on 
simulations with two degenerate dynamical 
flavours, which are identified with the u and d 
quarks. The strange and higher-mass quarks 
are still treated in the quenched approxima- 
tion. At this conference, CP-PACS reported 
that the resulting strange meson spectrum is 
much closer to experimented, as shown in Fig- 
ure ||. 

The glueball spectrum has only been 
computed in quenched QCD and the resultsE3 
reported at ICHEP98 remain state of the 
art. This calculation was hard because of 
strong scaling violations. Better scaling has 
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Figure 3. Comparison of continuum extrapolations 
of the K* and <j> meson masses in two-flavorp-i(filled 
symbols) and quenched (open symbols) QCDE3. The 
K mass was used to fix the strange quark mass. 



been re p orted recently using the fixed-point 
actionO. Mixing with quark states should be 
very important. A first attempt to compute 
the mixing of the lightest scalar glueball with 
the lightest scalar quarkonium states has con- 
cluded that the / (1710) is 74% glueball, 
whereas the /o(1500) is 98% quarkoniumEi 

Now that two-flavour simulations are 
possible, it is interesting to try to compute 
the flavour-singlet meson masses. A first at- 
tempt has produced a result that the 77, 77' 
mixing angle is around 45° in the {(uu + 
dd)/\ // 2,ss} basis, but, despite sophisticated 
variance reduction techniques, iit least ten- 
times better statistics is neededEa. 

5 Quark Masses 

Quark masses are encoded within hadron 
masses. They are scale and renormalisation- 
scheme dependent. To be useful for phe- 
nomenology, it is necessary to compute the 
scale evolution of the quark mass, from the 
lattice scale at which it is determined, to a 
suitably high scale where it can be matched 
to a perturbative scheme. 

The usual approach is to define an in- 
termediate scheme, such as the Schrodinger 
Functional (SF), in which the scale depen- 
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Figure 4. Continuum extrapolation of the RG- 
invariant mass, M a + M t , where M t = (M u + M d ) /2, 
in uuiis of the K decay constant for quenched 

qcdB 



dence can be com pu ted non-perturbatively 
(see Heitger's talkEZI). Once this scale de- 
pendence is known to a high enough energy, 
it can be continued to infinite energy, us- 
ing the perturbative renormalisation group, 
to define the ratio of the lattice mass to the 
renormalisation-group (RG) invariant mass, 
M§1. Thus the lattice quark mass fixes M 
and, from it, perturbation theory may be 
used to determine the quark mass in any cho- 
sen scheme. 

Figure [| shows the sum of the RG- 
invariant average u and d, and s quark 
masses, computed in this way fornuenched 
QCD (with the K mass as input The re- 
sulting estimate for the strange quark mass 
is 

mfS(2 GeV) = 97(4) MeV (N f = 0). (7) 

At this conference, CP-PACS announced 
results for light quark masses intwo-flavour 
QCEE3, updating earlier resultst£l. They use 
an improved action, with a fixed lattice size of 
(2.5 frn) 3 , and extrapolate downwards in the 
sea-quark mass, from masses corresponding 
to pseudoscalar-to-vector meson mass ratios 
above 0.6, to the average u and d quark mass. 
Currently, non-perturbative matching is not 
available for two-flavour QCD, so CP-PACS 
uses mean-field improved 1-loop matching. 
The results are shown in Figure |L for three 
different definitions of quark mass at non-zero 
lattice spacing, and compared with quenched 
QCD. The different definitions permit consis- 
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tent continuum extrapolations and the final 
results are from combined fits with a single 
limit. They find a big effect from the inclu- 
sion of dynamical quarks. The average u and 
d quark mass is reduced by roughly 25%: 

mjf (2GeV) = 3.44±°;gMeV(tf/= 2) (8) 

mjf (2GeV) = 4.36t^MeV(iV / = 0). (9) 

Treating the strange quark in the 
quenched approximation, CP-PACS finds 
that the 20% inconsistency in the strange 
quark mass in quenched QCD disappears 
with dynamical u and d quarks (within 10% 



errors), as can be seen in Figure 
tinuum estimates are 



The con- 



MS/ 



m~ J (2GeV) = 88IgMeV {K mass) (10) 

mf S (2GeV) = 90±^MeV (0 mass). (11) 

Again, the mass is reduced substantially com- 
pared to quenched QCD. Such a low strange 
quark mass suggests a large value of e'/e, 
and raises the interesting question how much 
lower the strange quark mass would be if 
it were treated dynamically. The result 
m s /m u d — 26(2)t3 agrees with the chiral per- 
turbation theory estimate of 24.4(1.5). 

The b quark mass, in this world in which 
only the u and d quarks are dynamical, ob- 
tained from the B s binding energy at leading 
order in 1/ ra^ and using NNLO perturbative 
matching, isEil 
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b 1 



4.26(9) GeV. (12) 
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Figure 6. Continuum extrapolation of the strange 
quark mass for Nf = 2 (filled symbols) and Nf = 
(open symbols) QCD, determined from the K 
mass jumpe r figure) and from the <j> mass (lower 
figure)!: 



6 Heavy-Quark Decays 

The calculation of hadronic matrix elements 
associated with the weak decays of b quarks is 



the most successful phenomenological aj 
cation of lattice QCD (see Kronfeld's talk 32 

6.1 Leptonic Decays and Mixing 



The top-quark CKM matrix elements and 
the neutral B q mass difference are re- 
lated through the hadronic matrix element 

fBa \ B B a ■ 



M^r,BSo{m 2 t /M^)M Bq 
x\V tq V t * b \ 2 flB Bg . (13) 



Am, = 



Traditionally, fs and Bb have been com- 
puted separately in lattice QCD. Quenched 
estimates for fg have stabilised in recent 
years, but suffer a relatively large irre- 
ducible scale uncertainty due to the quenched 
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Table 1. Summary of lattice results for B deq 
stants, presented by Hashimoto at Lattice 99 





N f = 2 


N f = 


Sb (MeV) 
fB. (MeV) 
Ib s 1 ]b 


210(30) 
245(30) 
1.16(4) 


170(20) 
195(20) 
1.15(4) 



approximation. Including two dynamical 
flavours increases fg by around 20%, al- 
though statistical errors currently overwhelm 
systematic effects. Table [l] shows the best es- 
timates of the B decay constants, presented 
at Lattice 99. The only direct comparison 
with experiment is for /rj s , and here the lat- 
tice estimates, 

f Ds = 241(30) MeV (N f = 0) 34 (14) 
f Ds = 275(20) MeV (Nf = 2)0, (15) 

are consistent with experiment (eg, ALEPH's 
result of 285(45) MeVEj), although the errors 
are too large to expose systematic effects. 

The combination Jb\/ Bb may be com- 
puted directly in lattice QCD and, to the ex- 
tent that systematic errors in fs and Bb are 
correlated, this may be more reliable than 
separate determinations of fs and Bb- A re- 
cent non-perturbathzEly-renormalised result 
in quenched QCD isE3 



/ b v/b b = 206(29) MeV. 



(16) 



Systematics should also cancel in ratios, so 
that quenched results such as23 

h_ 



0.74(5) 



h 



1.16(7) 



(17) 
(18) 



are probably the most reliable. 
6.2 Lifetimes 

Lifetime calculations are at an exploratory 
stage. Two groups have computed the 
B s lifetime difference, ATb b /^b s , obtaining 
O.O47(15)(16)0 and 0.107(26) (14) (17)11, to 



be compared with the experimental upper 
bound of 0.31. Although they use quite dif- 
ferent lattice techniques, the matrix element 
calculations are consistent, and the discrep- 
ancy in the final results is due to one using 
the quenched, and the other using the un- 
quenched value of fs, ■ 

The Ab lifetime is a puzzle, because 
the experimental measurement for the ratio 
r(A 6 )/r(B ) = 0.79(5) is significantly differ- 
ent from one, whereas, to leading order in 
the heavy-quark mass, all b hadrons have the 
same lifetime. A preliminary lattice calcu- 
lation, however, does indicate that spectator 
effects are significant at the 6-10% leveEJ. 

6.3 Exclusive Semileptonic Decays 

The main purpose of computing exclusive 
semileptonic B decay form factors is to ex- 
tract model-independent estimates of \V u b\ 
and \V c b\ from experiment. All the results 
I present were computed in quenched QCD, 
although dynamical-quaik effects are now be- 
ginning to be exploredEJ. They are expected 
to be around 10%. 

During the past two years, the most 
progress has been in calculating the form fac- 
tors for B — > ttIv, defined by 

M 2 - M 2 
(n(p')\V»\B(p)) = — 2 * q»f (q 2 ) 

+ (>+^- M| ~ M ' /) /+(g 2 )(19) 

q = p-p, V>* = jry fi v- (20) 

Here, lattice QCD fixes the normalisation of 
the form factors, unlike heavy-quark effec- 
tive theory (HQET). However, today's lattice 
spacings are too large to represent a high- 
momentum pion accurately. So the kine- 
matic range is restricted to near zero recoil, 
and model-independent results are only pos- 
siblc for the differential decay ratelfiJo, as 
shown in Figure (?]. Differential rates should 
be measured experimentally soon, at which 
point, direct comparison with the lattice re- 




Figure 7. Results for the B -j*-,tt£u differential de- 
cay rate,i-epmpiled by Lelloucfo, from FN A TrH and 
UKQCDcj. Vertical lines indicate the momentum 
range, 0.4 GeV< < 0.8 GeV, within which the 
lattice artefacts are minimised andr-epmparison with 
experiment should be most reliableq. 

suits will provide a model-independent esti- 
mate for \V u b\- 

Model-dependent extrapolation is 
needed to obtain the full kinematic range. A 
dipole fit to the UKQCD results for f+(q 2 ) is 
shown in Figure || (a simultaneous pole fit to 
/o(q 2 ) imposes the constraint /o(0) = /+(0)). 
This fit gives a total decay rate 

r/|Kfe| 2 = 9±l±i PS" 1 . (21) 

There is good agreement with the results 
from other groups and with light-cone sum 
rules at low q 2 . The form factor /o provides 
an important consistency check on the lattice 
results through the soft-pion theorem: 

MdaJ = fB/U (22) 

in the limit of zero pion mass. Whether this is 
satisfied by curre nt simulations is somewhat 
controversialljO, but this will have to be re- 
solved if we are to have full confidence in the 
lattice results. 

There have been no new results for B — > 
ptv and B — > K*j and the present status is 
described in Lellouch's reviewed. 

The form factors for the heavy-to-heavy 
decays, B — ► D^lv, are better suited to lat- 
tice calculations, because the recoil is smaller 



Figure 8. Results from various lattice groups for ihp 
B — » n£u form factor /_|_, compiled by LellouchEJ. 
The solid line is a fit to the UKQCD data and the 
dashed lines are light-cone sum rule results. 

and present-day lattices can cover the full 
kinematic range. However, HQET is able to 
determine the normalisation at zero recoil in 
the heavy-quark symmetry limit. So, lattice 
QCD is left with the tougher task of quantify- 
ing the deviations from the symmetry limit at 
physical quark masses, needed to extract \V c b\ 
from experiment, which requires few percent 
accuracy. The FNAL group has devised a 
technique for determining these power cor- 
rections at zero recoil, from ratios of matrix 
elements, in which many statistical and sys- 
tematic errors cancelE3cj, obtaining, eg, 

F D * (1) = 0.935(22)(f 1 )(8)(20). (23) 

From the recoil dependence of the form 
factors, it is possible to extract the Isgur- 
Wise function, £(wA^ UKQCD has done this 
in quenched QCDc3, using the B — » Dtv 
form factor. The resulting estimate of 
shown in Figure ^|, is independent of the 
heavy-quark masses, for masses around that 
of the charm quark, and is insensitive to the 
lattice spacings used. Thus, it is demon- 
strably the Isgur-Wise function for quenched 
QCD. 

7 Kaon Physics 
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FiguBe-i 9. The Isgur-Wise function for quenched 
QCEO, as a function of the recoil w. p is the slope 
parameter at lu = 1. 



7. 1 Mixing 

The B-parameter for neutral kaon mixing, 



B, 



3 (K \s 7f2 (l - l5 )ds^(l - l5 )d\K l 
8 <i? |s7M75d|0)<0| S - 7M 75d|^ ) 



(24) 

is probably the best determined weak matrix 
element in quenched QCD. The most reliable 
result comes from staggered quarks, because 
the mixing due to chiral symmetry breaking 
is non-trivial for Wilson quarks. The error in 
the continuum resulted, 



B K {2 GeV) = 0.628(42), 



(25) 



is mostly from perturbative matching, and 
should be reduced when non-perturbative 
renormalisation becomes available. Dynami- 
cal quark effects raise Bk by around 5% at 
fixed lattice spacing, but it is not known yet 
how this affects the continuum result. 

7.2 Non-Leptonic Decays 

Although there are no new results, there is 
renewed optimism for lattice calculations of 
K — > 7T7T decays (see Testa's talkS), be- 
cause we now have more sophisticated tech- 
niques, which may afford control over the 
severe cancellations between the matrix ele- 
ments concerned. The new chirally symmet- 
ric lattice formulations should avoid the mix- 



ing between operators of different chiralities, 
and the large measured value for e'/e reas- 
sures us that a signal should exist. 

The fundamental problem for lattice 
QCD is that, according to the Maiani- Testa 
no-go theorem, there is no general method 
for dealing with multi-hadron final states in 
Euclidean space. The traditional approach 
around this is to use chiral perturbation the- 
ory to relate those matrix elements which can 
be computed, such as K — > vacuum, n, or 
7T7T at unphysical momenta, to the desired 
physical matrix elementsca. A new proposal 
is to tune the lattice volume so that one of 
the (discrete) energy levels of the two pions 
equals the K mass, and then relate the tran- 
sition matrbcelement to the decay rate in in- 
finite volume£3. 

The impending flood of experimental 
data for non-leptonic B decays presents a 
formidable challenge to lattice QCD. Chiral 
perturbation theory no longer helps. Per- 
haps, the B — > 7T7T factorisationEil, proved for 
TOb ^ Aqcd, can be exploited in some way? 

8 Structure Functions 

Lattice QCD can provide the normalisation 
for parton densities. Ultimately, this should 
test QCD and the validity of perturbation 
theory. It enables us to disentangle power 
corrections and, where experimental informa- 
tion is scarce, such as for the gluon distri- 
bution for x > 0.4, lattice QCD can help 
phenomenology Dynamical quark effects are 
presumably crucial. Although results so far 
are for quenched QCD, this will soon change 
(see Jansen's talkcH). 

The traditional approach uses the oper- 
ator product expansion to relate moments of 
structure functions to hadronic matrix ele- 
ments of local operators: 



M n (q 2 ) = / dxx n - 2 F 2 (x,q 2 ) 
Jo 



= c(%7/^s(m))4 2) (m) 
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cf [Gey] 

Figure 10. The leading contribution, ni2, and the 
power correction, C2/g 2 , to the lowest non-trivial mo- 
ment of the unpola«ised nucleon structure function 
in quenched QCDEJ. The hatched strip is the re- 
sult without the inclusion of higher-twist effects (the 
width of the strip indicates the error). 



0(l/q 2 



(26) 



(2) 

The Wilson coefficients, C n , are determined 
in perturbation theory and the hadronic ma- 

(2) 

trix elements, A n are determined on the lat- 
tice. Renormalisation is the major source of 
systematic error, since the product C({i)A(p) 
must be independent of the scale /i. This 
is achieved using a non-perturbative interme- 
diate scheme, in the same way as for quark 
masses: 



A(»)=Z^)Z^(»a)A 



latt 



(a). (27) 



The INT=SF schemed uses a step scaling 
function to relate the matrix element, renor- 
malised at the lattice scale, to a high scale 
where perturbation theory can be used to de- 
termine the RG-invariant matrix element in 
the limit /i — > oo. This, in turn, may be 
related via perturbation theory to MS. At 
this conference, JansenEj reported that the 
average momentum of partons in the pion in 
quenched QCD, computed in this way, is 



(x)(2.4 GeV) = 0.30(3), 



(28) 



to be compared with the experimental result 
of 0.23(2), and confirming early lattice results 
that the quenched estimate is larger than ex- 
periment. 



A quite different method involves com- 
puting the current-current matrix element, 
(h\ J fx J v \h), which appears in the cross- 
section, directly on the latticeEj. The 
Wilson coefficients are determined non- 
perturbatively from matrix elements between 
quark states, by inverting 

(p\MQ)M-q)\p) =E^™L.^>>9) 

x(p|Ofc ) .. M „(a)b),(29) 

thereby avoiding mixing and renormalon am- 
biguities. Using 62 operators and 70 mo- 
menta to extract the C"s, and reconstructing 
(N\ J^J V | N) from them and nucleon matrix 
elements, QCDSF obtained the lowest non- 
trivial moment of the unpolarised structure 
functionEil, shown in Figure [l^. This indi- 
cates large power corrections and strong mix- 
ing between twist-2 and twist-4 operators. 

9 Machines and Prospects 

Progress in lattice QCD, and particularly its 
application to phenomenology, continues to 
be critically dependent on increasing com- 
puter power. Three machines dominate lat- 
tice QCD today. Historically, the first was 
CP-PACS's 300 Gflops (sustained) Hitachi 
SR2201. This has been operating since 1996 
and cost approximately $70/Mflops. The 
second was the QCDSP, custom built us- 
ing 32-bit digital signal processors, which has 
been sustaining 120 Gflops and 180 Gflops 
at Columbia and Brookhaven, respectively, 
since 1998. Its cost was around $10/Mflops. 
This year, APE's latest fully-customised 32- 
bit machine, called APEmille, began opera- 
tion at Pisa, Rome and Zcuthcn, sustaining 
around 70 Gflops in the largest configuration 
so far (this will double by the end of 2000). 
Its cost is $5/Mflops. These machines show 
an encouraging trend towards greater cost- 
effectiveness. 

In December 1999, an ECFA Working 
Panel concludedH that "the future research 
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programme using lattice simulations is a very 
rich one, investigating problems of central im- 
portance to the development of our under- 
standing of particle physics" . It also con- 
cluded that "to remain competitive, the com- 
munity will require access to a number of 10 
Tflops machines by 2003" and "it is unlikely 
to be able to procure a 10 Tflops machine 
commercially at a reasonable price by 2003" . 

Two new projects are targeting 10 Tflops 
64-bit machines, with a price/performance of 
$l/Mflops, by 2003. The QCDOC project, 
involving Columbia and UKQCD will em- 
ploy PowerPC nodes in a 4-dimensional mesh 
interconnect. The apeNEXT project, in- 
volving INFN and DESY, will continue the 
APE architecture of custom nodes in a 3- 
dimensional mesh. Two US projects, Cornell- 
Fermilab-MILC and JLAB-MIT, are explor- 
ing Alpha and Pentium clusters using com- 
modity (Myrinet) interconnect, in the hope 
that these commodity components can be 
made to scale to many thousands of proces- 
sors and that the intense market competition 
will drive the price very low. These develop- 
ments, together with the highly parallel algo- 
rithms employed for QCD, suggest there will 
be no obstacle to multi-teraflops machines for 
QCD except money! 

However, we still do not understand the 
scaling of our algorithms well enough to 
predict how much computer power will be 
needed. Our best estimates are that to 
achieve comparable precision to quenched 
QCD, in simulations with two dynamical 
flavours with masses around 15 MeV, will 
require between 15 and 150 Tflops years^J. 
However, we know nothing about simulations 
with light enough quarks for p — ► ttttI We still 
have a great deal to learn. 

In conclusion, the range of phenomeno- 
logical applications of lattice QCD continues 
to expand. Key developments have been im- 
proved actions (reported at ICHEP98) and 
non-perturbative renormalisation, both of 
which have considerably increased our con- 



fidence in matrix element calculations. Lat- 
tice QCD continues to drive the development 
of cost-effective high-performance computing 
technology and there is no technological limit 
in sight. The primary objective will be to ex- 
tend the range of quark masses which may be 
simulated reliably, and it is hard to see how 
we will get away with less than 100 Tflops 
machines. Finally, in the discovery of lattice 
chiral symmetry, we are witnessing the "Sec- 
ond Lattice Field Theory Revolution", and 
this will vastly increase the reach of ab initio 
computer simulations. 
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